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In near-infrared NaCo observations of the young brown dwarf 2MASS J0041353-5621 12, we discovered a companion a little less 
than a magnitude fainter than the primary. The binary candidate has a separation of 143 mas, the spectral types are M6.5 and M9.0 
for the two components. Colors and flux ratios are consistent with the components being located at the same distance minimizing the 
probability of the secondary being a background object. The brown dwarf is known to show Li absorption constraining the age to less 
than ~ 200 Myr, and it was suspected to show ongoing accretion, indicating an age as low as ~ lOMyr. We estimate distance and 
orbital parameters of the binary as a function of age. For an age of 10 Myr, the distance to the system is 50 pc, the orbital period is 
126 yr, and the masses of the components are ~ 30 and ~ 15 Mjup. The binary brown dwarf fills a so far unoccupied region in the 
parameters mass and age; it is a valuable new benchmark object for brown dwarf atmospheric and evolutionary models. 
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1. Introduction 
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Very low mass binaries are of particular interest for a number 
of reasons. First, the binary fraction at very low masses as well 
as their orbital properties carry important information about the 
way binaries form (e.g. Close et al. 2003; Burgasser et al. 2007). 
Second, all components of a multiple system share the same evo- 
lutionary history so that a comparison between binary compo- 
nents is free of a number of degeneracies. Finally, binaries offer 
a model-independent way to determine the mass through mea- 
surement of their orbital period. This third point is most impor- 
tant for our understanding of the evolution of low-mass stars and 
brown dwarfs in particular at young ages where models still have 
relatively large uncertainties. 

The number of very low mass binaries has grown rapidly 
over the last years. An updated list of binaries with total masses 
below 0.2 Mq can be found at www. vlmbinaries . org, this list 
carries 99 entries as of Jan 2010. Of particular importance are the 
brown dwarf binaries with independent age constraints because 
they deliver empirical constraints on brown dwarf evolution 
models. Furthermore, in order to determine the mass from or- 
bital motion on reasonable timescales, the orbital period should 
be short enough, and the binaries should not be too far away so 
that spectroscopic investigation is possible. Usually, young bi- 
naries are members of star forming regions that are located at a 
distance of lOOpc or more, which makes a detailed investigation 
of low-mass members very difficult. Therefore, young, nearby, 
low-mass systems are of very great value for our understanding 
of low-mass star and brown dwarf evolution. 

In this paper, we present the discovery of a new very low 
mass binary, 2MASS 0041353-5621 12 (hereafter 2M0041), that 
is nearby and probably very young. The age of 2M0041 is con- 
strained to be lower than ~ 200 Myr by the detection of Li 
in an optical spectrum CReiners & Basrii2009i) . iReinersi (l2009l) 
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presents evidence for accretion deduced from the intensity and 
shape of emission lines, in particular Ha. Ongoing accretion 
would suggest that 2M0041 may even be as young as 10 Myr. 
Space motion of 2M0041 is consistent with it being a member 
of the ~ 20 Myr old Tuc-Hor association or the ~ 12 Myr old 
y6 Pic association. Unfortunately, the distance to 2M0041 is not 
yet known because no parallax measurement is available so far. 
The distance from spectrophotometry would be 17 pc if the ob- 
ject was an old, single field star, but the real distance is of course 
larger because the young object is more luminous. As long as 
no parallax is measured for 2M0041, age, mass, and distance are 
free parameters that can be constrained by measuring the orbital 
period of the system. 



2. Data and Analysis 

Data were obtained with NaCo, the Nasmyth Adaptive Optics 
System (NAOS) and Near-Infrared Imager and Spectrograph 
(CONICA) at ESO's Very Large Telescope dLenzenet al.ll2003l; 
Rousset et al. 2003 ). Four images were obtained in service mode 
on August 14, 2009: One image was taken in 7, one in H, and 
two images were obtained in K^. We obtained ten individual 
frames, each of them being the average of 3 observations with an 
individual exposure time of 30 s each. Thus, the total exposure 
time is 15 min per image. Individual images were jittered using 
a 5" jitter box, this allowed efficient reduction of the sky back- 
ground. The adaptive optics system was used with the N90C10 
dichroic using 90% of the light for AO and only 10% for the sci- 
ence camera. We chose this option because no bright AO source 
was available nearby. All observations were carried out at low 
airmass (< 1.3) and good seeing conditions (< 0.8")- We used 
the S13 camera with afield of view of 14"xl4" and a pixel scale 
of 13.2mas/pix. For our analysis, we use the standard pipeline 
products that are provided for service mode observations. Data 
reduction includes dark subtraction, flat fielding, and sky sub- 
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Fig. 1. NACO K-band image of 2M0041, North is up and East is 
left. Axis labels denote coordinates in mas relative to the center 
of the image. 

traction using the jittered images. A close-up of one of the two 
/Tj-band images is shown in Fig.[T] 2M0041 is clearly resolved 
as a binary with two components that are somewhat different in 
brightness; the component to the SW appears brighter than the 
NE component. 

To measure the separation and the flux ratio of the two com- 
ponents, we estimate the PSF from the science images in an it- 
erative procedure. We start with an image of a standard star. We 
use the so-called zero-point images that are routinely taken for 
NaCo service observations. With this first guess of the PSF, we 
search for stars in our science images and determine their posi- 
tion and flux using the public domain IDL package StarFinder 
(Diolaiti et al. 2000). Next, we extract the PSF from the science 
image given the information on flux and position using the rou- 
tine ps£_extract from the StarFinder package. We use this 
PSF in the next step to iteratively search for the two components 
again and redetermine the shape of the PSF. 

The deduction of the PSF from the science image has the ad- 
vantage that we do not have to rely on a PSF that was taken at 
a diff'erent time, direction, and object brightness. The shape of 
the PSF sensitively depends on the adaptive optics performance 
which depends on the brightness of the reference target. On the 
other hand, the separation of the two components that we wish to 
distinguish is not much larger than the width of the PSF, partic- 
ularly in the 7-band. During determination of the PSF, this can 
lead to the problem that the algorithm constructs a PSF that con- 
sists of both components. To avoid this problem, we assume an 
axisymmetric shape of the PSF. After each PSF determination, 
we construct a rotationally averaged version of the PSF that we 
use for the next iteration of StarFinder. 

We successfully identified the two components in all four 
images using the procedure described above. As results, we ob- 
tain the positions and individual flux of the two components, and 
the PSF of each image. 

We show the separation of the two components, their flux ra- 
tio, and the width of the PSF from our four images in Table[T] 
The results for the separation are consistent with each other. The 
flux ratios from the two Ks-bwd images agree very well, and 
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Table 1. Separation and flux measurements in the four images. 



# 


Band 


Separation 


PA 


Flux ratio 


PSF width 






[mas] 






[mas] 


1 


J 


142.3 


44.90° 


2.19±0.11 


101.1 


2 


H 


142.4 


43.34° 


2.06 ± 0.06 


72.3 


3 


K, 


143.0 


43.46° 


1.88 ± 0.05 


48.2 


4 


Ks 


143.4 


43.25° 


1.88 ±0.05 


52.9 



the flux ratio seems to be a function of color, which is consis- 
tent with two components of different temperature. The width 
of the PSF also varies with wavelength, this is expected because 
the Strehl ratio is lower at shorter wavelengths (Strehl ratios are 
0.50 in Ks, 0.33 in H, and 0.18 in J). Note that in the K,,-ha.nA the 
width of the PSF is on the order of the diffraction limit. Going to 
even longer wavelengths would not lead to better image quality 
because the diffraction limit grows with longer wavelengths (for 
example, it is ~ 100 mas at L'). In the /-band, the two compo- 
nents are not very well separated. We tried to use a PSF con- 
sisting of different elliptical components, but did not succeed to 
produce results better than our procedure described above. The 
y-band position angle differs from the H- and /Tj-band results by 
about 1.5 degree, which can be attributed to the slightly irregular 
shape of the /-band PSF. In particular, irregularities in the wings 
of the 7-band PSF may lead to a small shift of the suspected 
stars' positions (here, 50.2 pix or ~2.5% of the 7-band FWHM). 
This probably leads to a slight dependence of the photocenter on 
the stars intensity in the /-band. 

To quantify the uncertainty of our measurements, we per- 
formed Monte-Carlo simulations in all three bands. In each 
band, we constructed a series of 1000 images of a binary con- 
sisting of two objects with different flux ratios, separations, and 
position angles. We uniformly varied these values around those 
found in the real data, and we used asymmetric PSFs consisting 
of two elliptical Gaussian components that scatter around a de- 
scription optimized to fit the rotationally symmetric PSF from 
the procedure described above. Noise was added to the artificial 
data in order to match the quality of the real images. In the H- 
and A'-bands, the fitting process turned out to be very stable. The 
clear separation of the two components due to the narrow PSF 
allows a robust fitting process even in the presence of an asym- 
metric PSF. We derive the uncertainties in the flux ratio from the 
scatter around the mean of our simulations. The 2cr-scatter of 
the flux ratio is 3% in H and 2% in Ks- The /-band uncertain- 
ties are larger The main reason is not only the larger FWHM 
of the PSF, but also the asymmetry in the PSF mentioned above. 
Thus, the fitting process in the 7-band image sensitively depends 
on assumptions on the shape of the PSF. Because deconvolving 
that shape is less reliable in the 7-band image, we decided not to 
simply adopt the /-band uncertainties from our Monte-Carlo ap- 
proach, but conservatively estimate them to be roughly a factor 
of two larger because of PSF-dependent systematics. We adopt a 
y-band flux ratio uncertainty of 5%, which well covers the scat- 
ter we found during our attempts using different PSF assump- 
tions. 



3. Photometry of the components 

Photometric measurements of 2M0041 in J, H, and Kg are 
avail able froni the Two Micron All Sky Survey (2MASS, 
Skr utskie et a"l]|2006h . The 2MASS magnitude of 2M0041 re- 
flects the combined flux from both components. With the NaCo 
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Table 2. Measured parameters of the system 





2MASS 


Individual photometry 


Parameter 


A+B 


A 


B 


J 


11.96 ±0.02 


12.37 ±0.03 


13.22 ±0.04 


H 


11.32 ±0.02 


11.75 ±0.02 


12.53 ±0.03 


Ks 


10.86 ± 0.03 


11.32 ±0.03 


12.01 ±0.03 




i.iU ± 0.04 


1.05 ±0.04 


1.21 ±0.05 


AJ 




0.85 ± 0.04 


AH 




0.78 


±0.02 


AK, 




0.69 


±0.02 


SpT 




M6.5 ± 1 


M9.0 ± 1 


Position Angle 




43.6° 


±0.6° 


Separation [mas] 




142.S 


±0.5 



observations, we can now separate the combined flux into the 
two components using the 2MASS magnitudes and the flux ra- 
tios f[ I fi in the three filters; for the 7-band we can write 



Ji = 7 + 2.5 log 



1 + 



(1) 



The results for individual J, H, and Ks magnitudes are sum- 
marized in Table|2] Under the assumption that the two compo- 
nents are indeed forming a binary (i.e., both are located at the 
same distance d), we can determine the individual spectral types 
from the magnitude differences and from the information about 
the combined spectral type from integrated light. Relations be- 
tw een absolute 7-magnitud e and Spect r al Typ es were reported 
by lDahn etal.1 ^02) and ICruz et al.l (l2003h . The combined 
spectral type of 2 M0041 is M7.5 with an uncertainty of +0.5 
(iPhan Bao & Bessel 2006). Using the linear relation from Dahn 
et al. we find that the Spectral Type difference between the two 
individual components is 2.5 spectral classes. The dispersion of 
objects around the relation from Dahn et al. is 0.25 mag, which 
translates into an individual uncertainty of 0.7 spectral classes 
for individual objects, and an uncertainty of 1 .0 spectral classes 
for the diff'erence between two objects. This dispersion is much 
larger than the uncertainty from the photometric differences of 
2M0041 A and B so that we can neglect the latter. 

The linear relation between Spectral Type and absolute mag- 
nitude also allows us to use the combined spectral type (M7.5) 
to anchor the spectral type range of 2M0041 A and B. Using 
flux ratio and individual magnitudes, we can construct an ar- 
tifical average magnitude of the combined system 2M0041AB, 
(Mj^AB - 12.63, M^ AB = 11.55), which in this system must 
coincide with Spectral Type M7.5 + 0.5. From that, we derive 
individual Spectral Types for 2M0041 A and B of M6.5 ± 1.0 
and M9.0 + 1.0, respectively. 

An independent check of our assumption that both compo- 
nents are indeed located at the same distance comes from the col- 
ors of 2M0041 A and B. J-K colors of late-M o bject s are given, 
e.g., inlHawlev et al.l (l2002h . iDahn et al.l (l2002h . and IWest et all 
( 120081) . The colors of 2M0041 A indicate a spectral type in the 
range M5-M9 while 2M0041 B falls in the range M8.5-early 
L. Thus, the colors and flux ratios of 2M0041 A and B support 
the assumption that both components belong to a binary system: 
Given a dispersion of 0.2 mag in the J-K relation and much 
smaller uncertainties in our measurements of flux ratios and col- 
ors, the difference in distance modulus between the two compo- 
nents is unlikely to be larger than 6 pc. 



Table|2] also contains the mean position angle and separa- 
tion. For the position angle, the individual measurements are 
weighted according to the inverse PSF widths of the images; the 
result is 43.6° + 0.6°. The mean separation calculated from the 
four exposures is 142.8 mas, we use the scatter of 0.5 mas as a 
conservative estimate for the uncertainty of this value. 

The proper motion of 2M0041 is ~ 140 mas yr"' 

(iPhan Bao & Bessel 20061), i.e., the primary travels roughly by 
one full observed separation per year in SB direction. The sec- 
ondary is still bright enough to be visible in archive nIR mea- 
surements if the separation to the primary was large enough. We 
can use nIR or red observations of the region to look whether 
any signatures of a second object can be found. Assuming that 
the secondary would show negligible proper motion, the objects 
should be separated by ~ 1.5" on the 2 MASS images taken in 
1999, which would be difficult to detect given the ~ 3" FWHM 
of the 2MASS PSF On the other hand, ESO.R-MAMA plates 
taken in 1988 should show two objects separated by about one 
FWHM, which also is 2.5-3". We found no second object at the 
position of 2M0041 in the 2MASS images from 1999 and in the 
ESO.R-MAMA images from 1988, and we see no signs of an 
elongated PSF which could be indicative of a barely resolved 
second object close to the primary (~ 1"). We conclude that it 
is very likely that the second object in our images is physically 
bound to the primary. Confirmation of binary status, however, 
can only be accomplished by verifying common proper motion 
in an exposure taken at a second epoch. 

4. System parameters 

From the new photometry, we can determine the parameters of 
the two components, and we can estimate mass and orbital pe- 
riod for a given age. The distance to 2M0041 can be estimated 
from the difference between absolute and apparent magnitude. 
For field objects, t he absolute magnit ude as a function of spe ctral 
type was given bv Dahn etal] (|2002|) and lCruz et"ai] (|2003|) . We 
can use that estimate for the lower Umit of the distance, which 
would apply if the objects were field stars. To calculate the dis- 
tance that would apply if the object was very young, we can then 
scale the absolute magnitude by the radius difference between 
young and old objects assuming that the temperature is not sig- 
nificantly changing with age fo r a given spectral type. W e use 
the radius-age dependence from Baraffe et al.l (Il998l 120021) . 

All parameters together with uncertainties are given in 
Table|3] Uncertainties are due to the spectral type uncertainty, 
which has an important effect on the absolute magnitude, and 
to this we add the uncertainty in J in quadrature. Under the as- 
sumption that 2M0041 A is a field star, the distance to the system 
is 24 pc. As expec ted, the new distanc e is larger than the 17 pc 
calculated earlier (iFahertv et al.ll2009l) . because in that calcula- 
tion, all the flux measured in J was assumed to come from a 
single (and cooler) object. We can now use the radius-age re- 
lation to estimate the distance for different ages of 2M0041 A 
(see .Reiners ,2009.) . The distances to 2M0041 for ages of 5 and 
lOMyr are 71 and 50pc, respectively. To estimate the semi- 
major axis, we apply a correction factor of 1.26 to account for 
projection effects (Fischer & Marcv 1992). The resulting sepa- 
rations between 2M0041 A and B come out as 8.9 and 12.8 AU 
for a lOMyr and a 5 Myr binary, respectively. 

From the evolutionary tracks of Baraffe et al., we can also es- 
timate the mass of the two components at a given age. With the 
masses and the separation known, we can then estimate the or- 
bital period of the binary assuming that the distance we observe 
is the true semimajor-axis of the system. 
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Table 3. System parameters for different ages. The main se- 
quence (MS) case is excluded by the Li detection and given only 
for reference. 



Parameter 




Value 




Age [Myr] 


5 


10 


MS 


dm 


71+-' 


JO_,(, 


24^^ 


Separation [AU] 


12.8!^-3 




4 3+1.8 


Mass (A) [Mj„p] 


-25 


~ 30 


-95 


Mass (B) [Mj„p] 


~ 15 


~ 15 


~ 80 


Orbital Period [yr] 




126!^« 






10 100 1000 

Age [Myr] 



Fig. 2. Top panel: Estimated orbital period as a function of age 
for 2M0041A+B (solid line). The grey region indicates Icr un- 
certainties (see text). Bottom panel: Estimated mass as a func- 
tion of age. The solid line shows the mass of component A, the 
dashed line is for component B, uncertainties are shown for com- 
ponent A only. The hatched area marks the region excluded from 
the presence of Li in the spectrum. 

The orbital period occuring if the binary was old is f ^ 22 yr. 
This is the lower boundary for the period, the real period must be 
longer because an age above a few hundred Myr is excluded by 
the detection of Li. If 2M0041 has an age of 10 Myr, the period 
on the order of 126 yr, i.e., a factor of 5 longer If the system is 
as young as 5 Myr, the period is about 228 yr, i.e., another factor 
of two longer. We show the estimated period as a function of age 
in Fig.|2] With the photometric information gathered about this 
object, the possible orbital period is a steep function of age. 

5. Summary 

In a NaCo image of the nearby young brown dwarf 2M0041, we 
have discovered the binarity of this object. The system consists 
of an M6.5 primary and a secondary of spectral type M9.0. We 



find a separation of 143 mas and derive individual J, H, and Ks 
magnitudes for both components. Flux ratio and colors of both 
components are consistent, which means that the chances for the 
secondary being a background object are very small. 

Distance, age, mass, and orbital period are yet unknown, but 
we can present possible solutions as a function of age. The object 
still has lithium, which means that it is a brown dwarf younger 
than a few hundred Myr. Ha and other optical emission lines 
indicate that the object may be accreting so that its age may be 
as low as 10 Myr For this age, the period is predicted to be on 
the order of 125 yr at a separation of ~ 9 AU and a distance of 
50 pc. Age, semi-major axis, and distance are steep functions of 
the orbital period. So far, no parallax measurement is available 
so that the distance to the binary is unknown given that we do 
not certainly know its age. 

Objects of spectral type M9.0 are much fainter at optical 
wavelengths than M6.5 so that the secondary contributes negli- 
gible background continuum to the Ha and Li measurement. On 
the other hand, the Ha emission seen in the combined spectrum 
may come from either of the two components, or both, which 
would affect the measured accretion rate; if both objects are ac- 
creting, the accretion rate per star would be lower, but the shape 
and strength of Ha still would indicate accretion. Alternatively, 
if the accretion timescale is longer on the secondary, it may be 
the only accreting object in the system. In that case, the accretion 
rate would be higher than the one originally derived. An image 
taken at Ha or even spatially resolved spectroscopy can solve 
this issue. 

The discovery of binarity in this young brown dwarf opens 
the opportunity to directly determine the mass of two young 
nearby brown dwarfs that may be accreting. Without knowing 
the distance, the age of the system can be determined from the 
orbital period. In the next years, the first estimate of the orbital 
period will become available. If it is indeed as young as a few 
ten Million years, it is the first object in this mass/age regime for 
which a direct mass estimate will become available. With a mea- 
sured parallax, independent information on the distance putting 
strong constraints on the system's parameters will be given. This 
system will be an important benchmark for brown dwarf evolu- 
tionary models at young ages. 
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